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Thermodynamic consistency for nuclear matter calculations 
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Abstract. We investigate the relation between the binding energy and the Fermi energy and between 
different expressions for the pressure in cold nuclear matter. For a self-consistent calculation based on 
a $ derivable T— matrix approximation with off-shell propagators the thermodynamic relations are well 
satisfied unlike for a G— matrix or a T— matrix approach using quasi-particle propagators in the ladder 
diagrams. 
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1 Introduction 

Nuclear matter calculations are usually performed using 
Brueckner type resummation of ladder diagrams. Works 
using realistic interactions lead to reasonable results for 

' the saturation density and the binding energy at the sat- 
uration point. However in violation of the Hugenholz-Van 
Hove theorem the resulting Fermi energy Ep at the satura- 
tion point is usually very different from the binding energy 

■ per particle E/N. It is a manifestation of a general viola- 
tion of thermodynamic consistency by the G— matrix ap- 
proximation. The problem was discussed in the literature 
M0 and improvements due to rearrangement terms were 
invoked but without removing the discrepancy altogether. 
Improvement of the fulfillment of the Hugenholz-Van Hove 
property with respect to the G— matrix approximation 
is observed when using the quasi-particle T— matrix ap- 
proach, or correction from hole-hole lines ]l|,0|. 

On the other hand it is known that the exact theory || 
||,H should fulfill certain thermodynamical relations. The 
simplest one being the exact equality of the Fermi mo- 
menta for the free and the interacting theory. Another 
statement that we shall consider in the present work is 
the equivalence of two ways of calculating the pressure in 
a system at zero temperature : 



P = p 



d(E/N) 



dp 

p(E F - E/N) , 



(1) 
(2) 



where p is the nuclear matter density. From the above 
relation follows that at the saturation point where (E/N) 
has a minimum 

E F = E/N , (3) 



i.e. the Hugenholz-Van Hove property. These relations are 
satisfied by the exact theory and can also be satisfied in a 
perturbative calculation to a given order of the expansion 
parameter. 

Non-perturbative approximations schemes which are 
thermodynamically consistent are know ^ . Baym has shown 
that the condition of the thermodynamical consistency 
of an approximation can be related to the so called 
derivability. The self-energy is constructed as a functional 
derivative of a functional of dressed propagators G(k) 
and bare vertices 

The approximate functional <P is defined by a set of two- 
particle irreducible diagrams. # derivable approximations 
to the self-energy are also termed as conserving approx- 
imations since they lead to conservation laws in corre- 
sponding transport equations |(|. In particular different 
types of non-perturbative approximations can be iden- 
tified for the generating functional. Below we shall con- 
sider two of them the Hartree-Fock approximation and the 
T— matrix approximation (Fig. |l). Diagrams for the cor- 
responding self-energies obtained by taking a functional 
derivative are also shown in Fig. |l[ It must be stressed 
again that the propagators in the diagrams for <P are dressed 
self-consistently by the self-energy (^). For the Hartree- 
Fock approximation it means only a shift in the single- 
particle energies, but for the T— matrix approach one has 
to take into account the full spectral function for the prop- 
agators in the calculation of <P or S. Calculations involving 
off-shell propagators in the T— matrix ladder have been re- 
cently performed [^,|^,[n],[l]J both in the normal and in the 
superfluid phase. Below we shall restrict ourselves to zero 
temperature normal nuclear matter. 
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Fig. 1. Diagrams contributing to the generating functional <2> 
in the Hartree-Fock approximation a) and in the T— matrix 
approximation c). The corresponding diagrams for the self- 
energy are shown in parts b) and d) respectively. 



2 Approximations for the nuclear matter 
problem 

We shall compare different calculation of cold nuclear mat- 
ter with a model interaction. We choose a separable rank 
two parameterization of Mongan type [l3| in the S wave 
with softened repulsive core 



V a (k, P ) = Kg r a (k)g r a ( P ) - Kg a a (k)g a a ( P ) 

with g r ' a (p) = pi ^i - and 

A r = 29.6GeV 2 (3 r = 639MeV 

\ a = 2.91GW 2 (3 a = 352MeV for a = S£ 

A r = h.27GeV 2 (3 r = 47lMeV 

\ a = 4.78GeV 2 (i a = 376MeV for a = Sf 



(5) 



(6) 



With this interaction nuclear matter properties will be 
calculated within the following approximations 

— Brueckner resummation of particle-particle ladder di- 
agrams with in medium G— matrix 



< P \G(p,n)\p >=v(p, P ) 

+ 



d 3 Q y(p)q) (l-/K 1 ))(l-/K 2 )) 



(2tt) 3 
<q|G(P,tf)|p' > , 



Q 



(7) 



where px,2 = P/2 ± q. G— matrix resummation al- 
lows to define single particle energies and gives rel- 
atively good results for the saturation properties of 
nuclear matter. In the above equation and in the fol- 
lowing we skip the spin, isospin indices which are im- 
plicitly summed over. Medium effects enter through 
the Pauli blocking factors 1 — f(u p ) in the numerator 
and single-particle energies to p in the denominator. The 
single particle energies u) p , are self-consistently defined 
by the G-matrix 



P 



(8) 



where 



U(p,u) 



(2tt) 3 



< (p-k)/2|G(|p + k|,^.+^)|(p-k)/2 > .(9) 

In the quasi-particle T— matrix approximation [[l2| . [Llfl 
the ladder diagrams include both particle-particle and 
hole-hole propagation. The Pauli blocking factor (1 — 
f(u> pi ))(l — f{uj p2 )) in the G-matrix equation is re- 
placed by 1 — f(u) Pl ) — f{w P2 ) in the equation for the 
retarded T-matrix 

<p|T(P,fl)|p' >=V(p,p') + 

3 v l^p, CjJ 



(27r) 3 ' n-u pi - uj P2 

<q\T(P,f2)\p > . 



(10) 



The imaginary part of the retarded self-energy in the 
T-matrix approximation is 

ImZ(p,u>) = J ^(/(c fc ) + 6(c + Wfe )) 

< (p - k)/2|ImT(|p + k\,LJ k + lu)\(p - k)/2 >(U) 

where b(ui) is the Bose distribution. The real part of 
the self energy consists of the Hartree Fock self-energy 
and a dispersive contribution obtained from IraE 



f duj 

RcS( P ,uj) = S HF {p)+V I — 



duj ImU(p, 



(12) 



The imaginary part of the self-energy is neglected lead- 
ing to the quasi-particle approximation for the two- 
nucleon propagator in the T-matrix (Eq. |Io|) . 
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Allowing for off-shell propagation of nucleons and tak- 
ing the self-energy self-consistently (also its imaginary 
part) requires the use of full spectral functions in the 
calculation resulting in more complicated expressions 
for the T— matrix and the self-energy |8j|j| 



<p|T(P,fl)|p > 

du>i f dw2 , 
~2T J IhT J (2tt) 3 



V(p,p 

d 3 q 



(1 - /K> - /(wa)) 



n — uj\ 

AfawdAfau*) <q|T(P,/2)|p' > 
and 

du>i 

< (p-k)/2|ImT(p- 
f/(«i) + 6(w 



■ w 2 



(13) 



Imi7(p, w) 



(2^ A(fc ' Wl) 
•k,w + Wi)|(p-k)/2 > 



(14) 



^ImZ^p, tu) 



p 2 /2m-Re£( P ,u)) +Im£(p 



where 



A(p,oj) 



(15) 

is the self-consistent spectral function of the nucleon. 
— Finally we present results for a simple Hartree-Fock 
approximation. It is certainly not well suited for real- 
istic applications in nuclear matter. However this ap- 
proach is (p derivable and it is illustrative to check 
its thermodynamic consistency explicitly. The Hartree- 
Fock approximation with parameters given by Eq. (|J) 
shows no saturation. We reduced the repulsive part of 
the interaction \ r a by 1.15 for the Hartree-Fock calcu- 
lation. This rescaling mimics the effect of ladder re- 
summation which leads to a reduction of the repulsive 
core. 

Equations for all the approximations schemes have to be 
solved iteratively, with a constraint on the total density. 
The numerical method for the solution of the T— matrix 
equation with off-shell propagators |l(| has been general- 
ized to the case of low and zero temperature. The details 
of the numerical procedure will be given elsewhere. 



3 Results for thermodynamic properties 
around the saturation point 

Only within the self-consistent T— matrix calculation 
the momentum distribution of nucleons 



i(p) 



M did ., . 
^-A(p,to) 



(16) 



55 I 



j=E f 



.7 at the Fermi momentum. 
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Fig. 2. Momentum distribution of nucleons for the full 
T— matrix calculation compared to the free fermion distribu- 
tion. 
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Fig. 3. Fermi momentum as obtained in the full T— matrix 
calculation (points) compared to the Fermi momentum of the 
free fermion gas (solid line). 



In the calculation the chemical potential /i = Ep is fixed 
by the constraint on the total density 



different from the Fermi-Dirac distribution (Fig. |2j). Clearly 
a Fermi liquid behavior is observed in the T— matrix ap- 
proximation, with a jump in the fermion density of 
-l 



2^ 



d 3 p 



A(p,w) = p 



(17) 



The corresponding Fermi momentum pp is defined by 
Ep = lu Pf . For a conserving approximation the Fermi 
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Fig. 4. Binding energy per particle for the G— matrix calcula- 
tion (solid line), the on shell T— matrix calculation (dashed 
line), the full T— matrix calculation (dotted line) and the 
Hartree-Fock calculation (dashed-dotted line) as function of 
the density. The corresponding Fermi energies are denoted by 
the same lines with solid boxes, open boxes, full circles and 
stars for the G— matrix, on shell T— matrix, full T~ matrix and 
Hartree-Fock results. In the insert is show a blow up of the 
region around the saturation point for the full T— matrix cal- 
culation. 



momentum should be the same as the Fermi momentum 
of a free fermion gas ffyM. Indeed it is well satisfied for 
a range of densities for the self-consistent T— matrix cal- 
culation (Fig ft) . All the other approximation discussed in 
this work fulfill this relation trivially since they use quasi- 
particles. 

The energy per particle in the different approximations 
can be obtained from the standard form of the energy 
density 



E/N= 1 - f**-*!Ll(£. 
' p J (2tt) 3 2tt 2 \2m 



+ u)A(p,u)f(u>). (18) 



Only for the self-consistent T— matrix the spectral func- 
tion and the u> integration is nontrivial. For the other ap- 
proximation schemes the spectral function is a delta func- 
tion. In that case the energy per particle can be expressed 
in the usual way using the single particle potential and 
kinetic energies. 

In Fig. H is plotted the energy per particle for differ- 
ent approximations for a range of densities aground the 
saturation density. The G— matrix and the full T— matrix 
calculations give very similar results for the binding en- 
ergy. The T— matrix with quasi-particle propagators gives 
somewhat different results, with lower saturation density 
and smaller binding energy. This behavior is due to very 
strong modifications of the effective mass around the Fermi 
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(MeV) 


K~ 
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(MeV) 


Hartree-Fock 


1.55 


-3.5 


-3.5 


87 


G— matrix 


1.42 


-21.6 


-10.9 


107 


T— matrix on shell 


1.08 


-18.9 


-7.0 


103 


T— matrix 


1.39 


-9.9 


-9.9 


103 



Table 1. Saturation density, Fermi energy, binding energy and 
compression modulus for different approximations discussed in 
the text 



momentum in the quasi-particle T— matrix approach. This 
effect is caused by the appearance of the pairing singular- 



ity in the T— matrix 14, 15|,|10[. In fact the quasi-particle 
T— matrix approximation is oversensitive to the presence 
of the pairing singularity, since the use of full spectral func- 
tions reduces the influence of the Cooper pair bound state 
on the nucleon spectral function and the single particle 
energies |l0|,|l6|]. 

In the Table are shown the corresponding binding en- 
ergies and saturation densities. The Hartree-Fock approx- 
imation gives significantly different results, only after a 
change in the parameters it has a saturation point at all. 
The Fermi energy obtained for different densities depends 
very much on the approximation chosen. Only for con- 
sistent approaches, i.e. Hartree-Fock and self-consistent 
T— matrix is the Hugcnholz- Van Hove condition at the 
saturation point satisfied. The difference between Ep and 
E/N at the saturation point is zero within numerical accu- 
racy for the Hartree-Fock and the self-consistent T— matrix 
calculations, and becomes as large as 10.7MeV for the 
G— matrix approximation. As previously observed the use 
of on-shell (i.e. non self-consistent) T— matrix approxima- 
tion ||] instead of the G— matrix reduces the violation of 
the Hugcnholz- Van Hove theorem at the saturation point. 
However, in the case where the pairing effect is strong (as 
in this work) the use of the quasi-particle T— matrix does 
not cure the violation of the Hugcnholz- Van Hove property 
and moreover deforms the results for the binding energy 
and the effective mass. We observe a very good fulfill- 
ment of the Hugenholz-Van Hove condition in the numer- 
ical solutions for <P derivable approaches with self-energies 
sclf-consistently taken into account, which means for the 
T— matrix calculation the use of self-consistent spectral 
functions in the propagators. 

The pressure can be calculated for a range of densi- 
ties by two methods (Eqs. |l|, ||) which should be equiva- 
lent. However only for the consistent approximations we 
find an approximate equivalence between the two formu- 
las, with very good agreement for the Hartree-Fock cal- 
culation (Fig. H). On the other hand non-consistent ap- 
proaches give very different results. In particular the point 
where the pressure equals zero and the slope of the pres- 
sure versus density comes out differently for the two ways 
of calculating the pressure. The slope of the pressure as 
function of density defines the compression modulus of 
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P (fm" 3 ) 

Fig. 5. Pressure obtained as a derivative of the binding energy 
(Eq. ^) for the G— matrix calculation (solid line), the on shell 
T— matrix calculation (dashed line), the full T— matrix calcu- 
lation (dotted line) and the Hartree-Fock calculation (dashed- 
dotted line) as function of density. The corresponding pressures 
obtained using Eq. ^| are denoted by the same lines with solid 
boxes, open boxes, full circles and stars for the G— matrix, on 
shell T— matrix, full T— matrix and Hartree-Fock results. 



nuclear matter 



op 



(19) 



which should be positive at the saturation point as a con- 
dition of stability. 

As expected p"^ ] non-conserving approximations give 
reasonable results for the binding energy and not for the 
Fermi energy. Thus, one should use Eq. |l| for the calcu- 
lation of thermodynamic properties (In particular JC — 
9^(p 2 ^(^))). The compression modulus for different 
approximations is given in the Table. Its value is similar 
for different approximations using ladder resummation. 
The values of K, obtained are smaller than in usual nu- 
clear matter calculations because we reduced the strength 
of the repulsive core. 



4 Conclusion 



T— matrix approach, i.e. when using the T— matrix with 
on-shell quasi-particle propagators. The full T— matrix and 
the G— matrix calculations give similar results for E/N. 
The binding energy is a physical result that can be used 
for the calculation of the pressure or compression modulus 
also in the non-conserving G— matrix approach. The same 
is not true for the Fermi energy which for non-consistent 
approaches is unreliable and leads often to unphysical re- 
sults if used in the thermodynamical relations. We note 
that the use of the quasi-particle approximation in the 
T— matrix resummation can lead to wrong results for the 
the binding and Fermi energies if the effect of pairing is im- 
portant. The self-consistent T— matrix and the G— matrix 
calculations are less sensitive to the fact that we have ne- 
glected the superfluid transition for cold nuclear matter. 

This work was supported in part by the Polish Committee for 
Scientific Research (KBN) under Grant 2P03B02019. 
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We have investigate the thermodynamical consistency of 
different approximations for nuclear matter. The <P deriv- 
able T— matrix approximation with off-shell propagators 
is conserving and fulfills these relations. The same is true 
for the simple Hartree-Fock approximation. On the other 
hand the usual G— matrix approximation violates badly 
the Hugenholz- Van Hove relation for pressure at zero tem- 
perature. The disagreement is reduced but not cured com- 
pletely when using a simplified version of the conserving 



